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Abstract 

A simple monatomic system in two dimensions with a double-well 
interaction potential is investigated in a wide range of temperature by 
molecular dynamics simulation. The system is melted and equilibrated 
well above the melting temperature, and then it is quenched to a tem- 
perature 88% below the melting temperature Tm at several cooling rates 
to produce an amorphous state. Various thermodynamic quantities are 
measured as a function of temperature while the system is heated at a 
constant rate. The glass transiton is observed by a sudden increase of the 
energy and Tg is shown to be an increasing function of the cooling rate 
in the preparation process of the amorphous state. In a relatively-high 
temperature region, the system gradually transforms into crystals, and 
the time-temperature-transformation(TTT) curve shows a typical nose 
shape. It is found that the transformation time to a crystalline state is 
the shortest at a temperature 14 ~ 15% below the melting temperature 
Tm and that at sufficiently low temperatures the transformation time is 
much longer than the available CPU time. This indicates that a long-lived 
glassy state is realized. 



1 Introduction 

For a long time, a lot of studies have been made in understanding the glass tran- 
sition by experimental, theoretical and computational methods [TJ |5J [5] . One of 
the difficulties in understanding the glass transition arises from the fact that 
almost all glass forming materials consist of many constituents. Even in a 
computer simulation, more than two components are needed in order to avoid 
crystallization[31[S]. The problem in a simple monatomic system is that crystal- 
lization occurs in a relatively short time. Since a long time simulation is needed 
to investigate the slow dynamics around the glass transition point, a binary 
mixture is usually investigated. 

In order to avoid the crystallization of simple liquids in simulations, the 
Dzugutov potential was proposed [6^ . The Dzugutov potential has an additional 



maximum to the Lennard- Jones potential at a range typical of the next-nearest- 
neighbor coordination distance in close-packed crystals, which suppresses the 
crystallization. Indeed, the system keeps a metastable supercooled state for a 
long time and shows some characteristic behaviors of glass- forming liquids [21 E]. 
However, it turned out to form a dodecagonal quasicrystal by freezing in a 
subsequent simulation [5]. 

In fact, features of a glass formation are incompletely understood even for 
the simplest system, so it is attractive to carry out the study using a sophisti- 
cated monatomic system. In multicomponent systems, the relaxation involves 
both the topological and the chemical ordering. A spherically interacting one- 
component model system provides an opportunity for separating these contribu- 
tions and enables us to explore the topological effect to the relaxation dynamics 
around the glass transition point. It is useful for comparison with theories, since 
most theoretical works [9} fTO] ITT] assume a simple system. 

For the purpose of studying the glass formation of a simple liquid, we em- 
ploy a monatomic model system in two dimensions with a double-well interac- 
tion potential, the Lennard-Jones-Gauss(LJG) potential [TH IS]. In a previous 
report [2], we showed that the LJG system can be vitrified in two dimensions 
and the glassy state at low temperatures is stable for a fairly long time in spite 
of a simple monatomic potential. In addition, for 3D systems the glass-forming 
ability of the LJG system has been tested and discussed [15]. 

In this paper, we show the further study of the glassy feature in this LJG 
system, in addition, we investigate the aging-induced crystallization for the 
monatomic glass-forming system with MD simulations. In general simulations, 
the model system is chosen to favour the clarification of dynamical relaxations 
under supercooling, while the crystallization is highly suppressed. Therefore, 
correct discriptions of the relaxation process of the glass-forming liquids are 
still not really understood, and we focus attention on this point. In Sec. 2, we 
explain our model and methods of MD simulation in detail. We present the 
evidence of vitrification of this system in Sec. 3 and we show the results of the 
crystallization process in Sec. 4. Finally, we discuss the results in Sec. 5. 

2 Model Potential 

We consider a system of atoms which interact isotropically through the Lcnnard- 
Jones-Gauss (LJG) potential. 

The LJG potential consists of the LJ potential and a pocket represented by 
the Gaussian function. It is a double-well potential for most values of the 
parameters with the second well at position re, of depth e, and width a. We 
note that this kind of effective potential is known to exist in matals. The general 
form of pair potentials in metals consists of a strongly repulsive core plus a 
decaying oscillatory term [T6|. The LJG potential can be understood as such 
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Figure 1: The LJG potential for the parameters tq — 1.47ro,e — 1.5eo,cr^ = 
0.02r§ 

an oscillatory potential, cut off after the second minimum. This potential was 
proposed for the self-assembly of two-dimensional monatomic complex crystals 
and quasicrystals [T^]. It has also been used as a model of liquid water [17] and 
as a potential to stabilize a given structure [TH'. We find that stable glassy states 
can be obtained by quenching simulation with several values for the parameters. 

In this paper, we fix the parameter tq — 1.47ro, e — 1.5eo, = 0.02rQ. Fig- 
ure [l] shows the LJG potential with this chosen parameters. This set of param- 
eters produces the pentagon-triangle phase in the ground state |13]. The unit of 
length, energy and time in the present simulation are rg, eoi t = ('^trQ/eo)^^^, 
respectively. 

We perform MD simulations for 2,500 atoms on two dimensional space with 
periodic boundary conditions. The size of the simulation cell is chosen so that 
the number density becomes p — 1.0. We employ a Verlet algorithm with time 
step O.OIt. 

3 Results and Discussions 
3.1 Glass Transition 

In order to determine the glass transition point, we follow the procedure em- 
ployed in the thermodynamic measurement [H]. Namely we prepare an amor- 
phous state at a far below Tm with a rapid cooling, and anneal the system for 
some time. Then we heat the system at a constant rate and measure the energy 
of the system. The glass transition temperature is identified by a change in the 
slope of E vs T graph. 

We prepare an equilibrium liquid state as the initial condition at T* = 
ksT/eQ — 2.2 which is above = 0.43 and then cool the system down to a 
target temperature with the velocity scaling. When the target temperature is 
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Figure 2: Tiling structures in real space and diffraction images of a glass and 
crystal after a long run. A blue and yellow tile show pentagon and triangle, 
respectively. Square and hexagon are painted by same color, red. A gray tile 
shows collapsed nonagon. 

close but below T^, the system transforms into a crystalline state after some 
time. However, the target temperature is well below Tm, the system remains in 
the amorphous state even for a long time run. Figure [2] shows tiling structures 
in real space and diffraction images of a glass and crystal after 3.5 x lO^r. We 
determined tilings by drawing a line between the particles within a cutoff radius, 
where the cutoff radius is chosen as the first minimum of the pair correlation 
function of the system. Triangle and Pentagon tiles, which are painted by yellow 
and blue respectively, are essential for stable structures of this system. In the 
glassy state, any long range order cannot be observed. 

Figure |3]shows the radial distribution function of a glassy state. The posiiton 
of the first peak is 0.93 and the second peak is 1.50. The position of the second 
peak approximately coincides with the second minimum of LJG potential. The 
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position of the first peak corresponds to a slightly inside from the first minimum 
of the potential. Here, we consider a regular pentagon, and if the distance 
between second nearest particles has the value of 1.47, the distance between 
nearest particles becomes about 0.91. This means the potential used in this 
study favors a pentagonal local order especially because of the position of the 
second minimum. This pentagonal local order is essential for the stability of the 
glass in this system. The third and fourth peak represent the existence of the 
local structure consisting of one pentagon and one triangle(see Fig. [s]). These 
tiles share on edge and make an ordered arrangement in a crystalline state, 
while it randomly spreads out all over the system in a glassy state (see Fig. |2|. 

Figure |4] shows the change of the total energy in a heating process from a 
glassy state at T* = 0.05 prepared with cooling rate = 10~^eo/kBT. The 
heating rate from a glassy state is jh — lO^'^co/fcsT. The slope of the black 
line in Figure |4] is about 2, which represents the energy change due to the 
vibration. Simulation data lie on the black line at lower temperatures, which 
indicates that the specific heat is determined by the vibrational motion of atoms. 
However, plots start to deviate from the line at T* ~ 0.38, which is below the 
melting temperature, so that an additional effect appears in addition to that 
of vibration. Crystallization does not occur in this heating process, which is 
confirmed by the static structure factor. That means the system transforms 
from a glassy state into a supercooled state as the temperature is increased. If 
we take -fh = lO^'^eo/^s''" as the heating rate, the sign of crystallization appears 
in the heating process, so we cannot decide the glass transition point with 
confidence. Thus, we conclude the glass transition temperature is T* ^ 0.38 in 
this sample. We perform the heating simulation starting from glassy samples 
prepared by different cooling rates and determine Tg for each sample. Figure 
[5] shows the cooling rate dependence of Tg . We find that Tg is an increasing 
function of the cooling rate, which is consistent with experiments |20j . 
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Figure 4: The temperature dependence of the total energy in a heating process 
from a glassy state. If the system can be described only by the term of lattice 
vibration, plots are supposed to lie on the line in above figure. 




cooling rate [10 So/kgi] 

Figure 5: The cooling rate dependence of Tg 
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3.2 Time- Temperature- Transformation Diagram 

In this section, we investigate the structural transformation into a crystal in 
the real space. We first prepared an equilibrium liquid state at T* = 0.4 above 
the melting temperature ~ 0.34. Then, we quenched the system to zero 
temperature instantaneously by removing the kinetic energy. We employ the 
instantaneous cooling in order to avoid any nucleation during the cooling process 
which may hinder proper comparison of the crystallization process. After the 
quenching, we subjected the system to a heat reservoir at the target temperature 
and allowed it to relax at this temperature. The simulation cell is made large 
enough so that the system feels zero pressure at all times. Under this free 
boundary condition, particles on the surface can easily move, so rearrangement 
of particles can frequently occur, compared with periodic boundary conditions. 
Since crystallization proceeds more drastically, we can determine easily the time 
when the system crystallizes. The number of atoms in the system is = 1, 024 
and the Leap- Frog algorithm is employed. 

We checked the time needed for crystallization at each temperature and de- 
termined a time-temperature-transformation(TTT) diagram. In order to deter- 
mine the crystallization time, we focus on the formation of three characteristic 
tiling structures, PPPT, PPP and PPTS, which consist of pentagons, squares 
and triangles shown in Fig. |6|a). Figure [6j^b) shows the time dependence of the 
number of atoms whose surroundings can be identified as the three characteris- 
tic tiling structures. The curve (all) shows the number of particles surrounded 
by any of the pentagons, triangles or squares. This value is related to the poten- 
tial energy and becomes constant when the system reaches a crystalline state. 
At relatively high temperatures, the transition from the supercooled liquid to 
a crystal is rather sharp and we can determine the crystallization time without 
ambiguity. As the temperature is reduced, the transformation tends to be less 
sharp, and we determined the transition time by observing diffraction patterns 
as well as the evolution of tilings. 

Figure |7] shows the crystallization time as a function of temperature which 
is known as a time-temperature-transformation (TTT) diagram. This diagram 
indicates that the crystallization time becomes the shortest at ^ 0.86Tm and 
it is longer than lO^r at 0.5Tm- If one uses the value of m, rg and eo relevant 
for the LJ potential of Ar, this time, lO^r corresponds to 10 micro seconds. 
The glass transition temperature Tg shown in Figure |7] is determined under 
the following condition. First, we prepared glassy samples by instantaneous 
cooling to O.lbTm- After that, these samples were annealed at this temperature 
for lOr, lO^r, lO'^r, lO^r, lO'^r and we used these 5 samples as initial states for 
observation of glass transition. The temperature was increased from these initial 
states by 10~^eo/kB a step and the energy was determined by its average for It 
period at each temperature step. Finally, we determined Tg by a sudden increase 
in E vs T graph. We found that Tg does not depend on the waiting time and is 
located at the temperature which a crystallization time is the shortest. 

For comparison, the results for a monatomic Lennard-Jones(LJ) system are 
also plotted in Fig. [?) Consequently, the LJG system has a much longer 
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Figure 6: (a) Three charateristic tiling structures, (b) The time dependence 
of the number of characteristic tiUng structures. The curve (all) represents the 
number of atoms surrounded by any of the pentagon, triangle or square tiles. 
The curves (PPPT), (PPP) and (PPTS) correspond to the structures shown in 
(a). 

crystallization time than the LJ system, and we can clearly see the temper- 
ature dependence of the TTT diagram in the LJG system. It shows a typical 
nose shape, which has been observed by experiments for various glass forming 
materials IHJ [551 [23]. A similar TTT diagram is also found by MD simulations 
with some empirical potentials for metal[211llS]. Our system has a much longer 
crystallization time than that reported in these papers in spite of the simpler 
shape of our interatomic potential. 

A LJ liquid rapidly transforms into the triang le lattice at a few 10^ MD 
steps. At low temperatures, particles condense into some clusters which are 
separated by some voids, and then such voids are eliminated gradually. On the 
other hand, at high temperatures, the system can transform in such a way that 
voids do not appear because particles can move more freely. This is the reason 
that the crystallization time in a LJ system becomes slightly longer at higher 
temperatures. 

Figure [8] shows results for the relaxation time dependence of potential energy 
at temperatures from 0.338 to 0.32 and from 0.23 to 0.20. It can be found 
from Figure [sja) that the system at temperatures very close to Tm{= 0.34) 
behaves as if it is an equilibrium liquid state for a while, and it drastically 
transforms into a crystalline state by random nucleation. The crystallization 
time in this temperature region is distributed because it depends on whether 
nucleation occurs or not by chance. Figure [9]shows one example of the process of 
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Figure 7: The time-temperature-transformation diagram for monatomic LJ and 
LJG liquids 

nucleation and growth at temperatures close to . The nucleation occurs in the 
inner part(Fig. |9|l) and grows rapidly(Fig. ^p), however in a while, the order is 
destroyed from the outer side (Fig. |9];). Since particles, especially on the surface, 
can move around rather freely, rearrangement can easily occur. Therefore, the 
system can remove defects and eventually reach a perfect crystal(Fig. [9|i). In 
Figure |8]d, the potential energy at a low temperature shows a higher value for 
long time scales compared with those at higher temperatures. It indicates that 
the system at sufficiently low temperatures keeps a metastable state for a long 
time. 

Figure [To] shows the system-size dependence of TTT diagram. The number 
of particles in each system is N =1024 and 4096. In both systems, a longer 
time is needed for crystallization near T^, and the temperature which has the 
shortest crystallization time is almost same. Tm in a larger system is higher 
than that in a smaller system. This tendency agrees with the size dependence 
of Tm found in confined systems [261. We also find that the shortest time in a 
large system becomes longer than that of small systems. It is because, in a large 
system, a much longer time is needed for the spread of crystalUne phase to the 
whole system. However, at lower temperatures, the growth of the crystallization 
time is suppressed in a large system. The reason of this result is that the state 
is different between two systems at the time when the system is assumed to 
have crystallized. At this time, many samples in the =4096 system are a 
polycrystalline state which has more than two domains, on the other hand, the 
N =1024 system has less domains and sometimes one domain spreads to the 
whole system. If we determine the crystallization time when the system has one 
or two domains, a large system will take a longer time to crystallize. 

We also investigated the crystallization time under periodic boundary condi- 
tions(PBC). The boundary condition dependence of the TTT diagram is shown 
in Figure [TT] The number of particles is 1024 in both systems. It is difficult 
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Figure 8: Time dependence of potential energy at several temperatures, 
shows the region close to T„i, 0.32< T <0.338. (b) shows the region at 
temperatures, 0.20< T <0.23 
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Figure 9: The process of nucleation and growth at a temperature close to Tm- 
(a)The nucleation occurs in the center and (b) it grows rapidly. (c)A part 
of crystalline phase is desroyed and rearrangement occurs. (d)Eventually, the 
system becomes an almost perfect crystal. 
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Figure 10: The system-size dependence of the time-temperature-transformation 
diagram for monatomic LJG liquids. The number of particles in each system is 
1024 and 4096. The temperature range in N=4096 system is from 0.24 to 0.341. 
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Figure 11: The boundary condition dependence of the time-temperature- 
transformation diagram for monatomic LJG liquids 

to determine the crystaUization time clearly under PBC in the same way as 
under free boundary conditions(FBC). In both systems, first nucleation occurs 
at almost same time. The difference appears in a subsequent process. In a 
system under FBC, if the nucleus starts to grow, a crystal phase spread rapidly 
to the entire system. It is more diflicult for the nucleus to grow under PBC 
and a perfect crystal cannot be obtained unlike FBC. We determine the TTT 
diagram under PBC by drawing upon the potential energy. Since the index of 
crystallization is different, the crystallization time cannot be simply compared 
with one under FBC. The similarities between these two systems are that the 
crystallization time becomes longer at temperatures close to Tm and it is the 
shortest at a temperature 15% below Tm- The difference appears in a behavior 
at low temperatures. The growth of the crystallization time in PBC system is 
suppressed than that of FBC system. Also in this case, the state is different be- 
tween two systems at the time when the system is assumed to have crystallized. 
A crystalline state in PBC system is less perfect than that of FBC system. It 
is more difficult in PBC system to get a perfect crystal because of boundary 
constraint. 

4 Conclusions 

We presented strong evidences of vitrification of a simple monatomic liquid in 
two dimensions by an adequate choice of LJG parameters. We have shown 
that the glassy state can survive a fairly long time at low temperatures and 
the glass transition temperature is an increasing function of the cooling rate 
in the preparation process of the amorphous state. The Gaussian part of the 
LJG potential stabilizes the pentagonal configuration and packing of pentagons 
produces frustration in crystallization. Two competing length scales and the 
local structure made by them can have a great impact on the degree of disorder 
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in the system, which leads to suppress crystallization. It has also been shown 
that, in 3D systems [T5]. the glassy state can be formed if the LJG potential has 
appropriate parameters favoring the formation of an icosahedral local order. 
The effect of this frustration with packing is the origin of the stability of the 
glassy state. 

One of the features of this system is that the time needed for crystallization 
is sufficiently long. It enabled us to determine the TTT diagram of the simple 
system. It has a nose shape and the transformation time into crystal is the 
shortest at a temperature 14 ~ 15% below , which are independent of the 
system size and boundary conditions. Above the glass transition temperature, 
crystal nuclei rapidly spread to the whole system. Below Tg , the growth of crys- 
tal nuclei becomes slower as a temperature is decreased, because rearrangement 
of atoms occurs less frequent. 

This model will enable us to explore in greater detail the mechanism of 
nucleation and growth processes in the monatomic glass. The study in this 
direction will lead us to solve the problem about slow nucleation kinetics of 
glass forming liquids. Due to simplicity of a model, one can use the LJG system 
for further study of other phenomena related to the glass transition and for a 
detailed comparison with theories. 
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